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We consider the model giving rise to Witten's superconducting cosmic strings at finite fermion 
chemical potential and temperature. We demonstrate how various symmetries of the hamiltonian 
can be used to exactly compute the fermion electric current in the string background. We show that 
the current along the string is not sensitive to the profiles of the string fields, and at fixed chemical 
potential and temperature depends only on the string winding number, the total gauge flux through 
the vortex and, possibly, the fermion mass at infinity. 



PACS numbers: 



I. INTRODUCTION 

Ever since Witten's pioneering paper [lj, it has been known that cosmic strings can posses fermion zero modes 
concentrated in the string core. One remarkable feature of this system, is that an application of a constant electric 
field in the string direction induces an electric current along the string carried by the zero modes. This current 
will grow linearly with time, while the electric field is turned on and will persist even after the field is turned off. 
The string, thus, becomes superconducting. It must be noted, however, that the behaviour of the system is known 
precisely only for induced currents smaller than a certain critical current - once the current exceeds this critical value, 
the energies of the zero modes become larger than the fermion mass at infinity m, and it becomes possible for the 
charge carriers to move off the string, quenching the superconductivity. The question of build upof charge and current 
on the superconducting string in an external electric field has been analyzed extensively in 0, |2, d,|il,0iE|- 

In this paper, we investigate a very different mechanism for inducing a current on the string. Namely, we compute 
the current on a superconducting string in the presence of a non-zero fermion chemical potential \x and temperature T. 
It is a rather trivial exercise to calculate the current J along the string for fi < m, T <C m, when then the low energy 
dynamics of the fermion-string system are governed by an effective 1 + 1 dimensional theory of zero modes moving 
along the string. In this case it is straightforward to show that the current for each fermion species is J — , where 
e is the fermion charge and n is the winding number of the string. However, we make a much stronger statement: it 
is possible to calculate exactly the total electric current in the string direction for any value of the fermion chemical 
potential fi and temperature T. We shall show that the result is topological in nature, and is independent of the 
particular profiles of the background string fields. The result will depend crucially on whether the string is local (as 
considered by Witten) or global (as, for instance, in the case of axion strings). In particular, if the string is local, the 
naive prediction for J of the effective 1 + 1 dimensional theory, remains valid for any value of fi, T. 

The appearance of quantum numbers (particularly of fermion number) on topological defects is a very well-developed 
subject with known computational methodsQ, such as trace identities and adiabatic expansion. At zero chemical 
potential, the fermion charge induced on defects is usually a topological quantity and, frequently, can be evaluated 
exactly. However, at arbitrary finite chemical potential, the fermion charge induced is generally not topological 8], 
and difficult to compute exactly. In view of this, our result is particularly interesting, since we show that quantum 
numbers such as total current can remain topological and exactly calculable at arbitrary fermion chemical potential. 
Mathematically, our analysis can be easily generalized to a large class of Hamiltonians involving fermions in d + 1 
dimensions in the background of a d dimensional defect, which is uniform in the [d + l)'st direction. 

As far as we know, the problem of computation of electric current on a superconducting string in the background 
of an arbitrary fermion chemical potential and temperature has not been considered before, although some of the 
techniques we use have been previously discussed in conjunction to index theorem 9] for string zero modes, and charge 
induced on the string by an electric field 4, 5j. Although the present problem could be of some interest in application 
to cosmology, this paper was largely motivated by related problems in dense quark matter. It is well known that 
quark matter at large baryon chemical potential, which might be realized in the cores of neutron stars, breaks certain 
symmetries of QCD^(J, and may supports several kinds of strings [llj. Recently, the method of fictitious axial 
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anomalies has been used 01 to derive the following effective action for the interaction of an electromagnetic field A h 
with an axial string in dense quark matter: 



E 



2tt ~q 



AidF (1) 



Here the line integral is along the string, the index a runs over all species of quarks, and the fraction ^ denotes 
the flavor content of the condensate that supports the axial string. This implies that the axial string in dense QCD 
carries an electric current of J = " for each quark species. We wish to understand this phenomenological 
result microscopically. The method for computation of current on strings at finite chemical potential developed in|l2| 
is sensitive only to the pattern of symmetry breaking and, thus, would yield in application to the present model the 
same result, J = As we see below, in some cases this result remains correct for all /i, while in other cases it 

receives corrections of order 1. 

II. CURRENTS ON STRINGS 
A. The Model 

Consider the following model of a Dirac fermion ip coupled to a string: 

C = ^(fy - ieA» - ^Rrfty " + <^^)V> (2) 

Here and i? M are gauge fields and <j> is a complex scalar field. The model has the following classical gauge 
symmetries: 

U(l) : V -> e^V, A»^A (l + d^a, 0^0 (3) 
U(l) : t/j -> eW&rfPi/), R v ,^R fi + 0^6, <p -> (4) 

This model is exactly equivalent to Witten's model of superconducting cosmic strings with a particular choice of gauge 
charges 1 . By convention, we associate the vector field A^ with electromagnetism. We note that the above model 

suffers from gauge anomalies, which can be removed, for example, by adding another fermion ip to the model with the 
opposite R charge q = —q and the electric charge e, such that e 2 = e 2 . The Lagrangian for the ip fermion is then: 

t = ^7"(0/i - ieA» + |^7 5 )^ - ht{4> l A^- + (5) 

Notice that ip now couples to <f>* rather than <f>. 

We could also consider the situation when the C7 (1) symmetry is global, such that the gauge field R^ is absent, the 
Lagrangian @ is by itself anomaly free, and the addition of the ip field is unnecessary. In our calculations, we will 
recover this case by taking q = 0. 

We assume that the U(l) symmetry is spontaneously broken, the <fi field acquires a non-zero expectation value and 
strings of the (j) field are possible. We wish to consider the fermion ip in the background of an infinitely long static 
string uniform in the z direction. The string is characterized by a non-zero winding number n of the scalar field: 

dl a cP*d a cj> 

2ni |0| 2 1 ' 

where the integral is over a contour in the xy plane at infinity, and the absolute value of the scalar field \<p\ tends to 
some constant (f> as r — > oo in the xy plane. 



1 We could have easily considered the completely general version of Witten's model, however, to simplify the algebra slightly we concentrate 
on the above choice of gauge charges, which makes the field couple to the vector current and the field couple to the axial current. 
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If the U(l) symmetry is local, then in most models (such as the Abelian Higgs model), the condition that the string 
energy is finite, implies that D^4> = (cL — iqRn)(f> — ► fast enough as r — > oo in the xy plane. This in turn implies 
the quantization of the string flux: 

$ = £ / d 2 xe ab R ab = n (7) 

From here on a, b = 1,2 and R^ v is the usual field strength tensor. It must be noted that the condition Q is not 
present in the case of global strings, so we will throughout our calculations keep the flux $ arbitrary and at the end 
set $ = n for local strings and $ = for global strings. 

Our objective is to calculate the expectation value of the electromagnetic fermion current in the string direction, 

J 3 = e J d 2 a#7 3 V>) (8) 

at finite fermion chemical potential and temperature T. Note that if the U(l) symmetry is local, there is an 
additional contribution to the electromagnetic current from the tp fermions. This, however, can be obtained from the 
result for the ip fermions by setting e — ► e, q — > q = —q, (p — > <p* , which translates into $ — > — <&, n — > —n. 

B. The Spectrum 

Let's start by analyzing the spectrum of our fermions in the string background. The one-particle Hamiltonian is: 

H = -ia l (d t - ieA t - |i? l7 5 ) + fc 7 °(0*i±X- + ^—^-) (9) 

where a 1 = ^y°j l and z = 1,2, 3. For a static background string uniform in the third direction, Aj = 0, R3 = 0, and 
hence, 

H = -id 3 a 3 + H 1 ^ (10) 
ff x = - l a a (d a - l -^R a ^) + h 1 ^^ l -±^+4> 1 —^-) (11) 

Since all the fields are assumed uniform in the third direction, we can choose —id^ip = kip and work at fixed k. 2 In 
each k sector, 

H k = ka 3 + H 1 - (12) 

and the operator now acts solely in the transverse xy plane. At this point, we make our choice of the 7 matrices 
to be: 



The operator H 1 - then takes the form: 



11 1 = ( vr u ) , l 4 ) 



where, 



2? = a a a a + |i? ae a V + /i(i^0+^^0*) (15) 



2 We take the third direction z to be compact of length L so that the eigenvalues k are discrete. As usual, we will take L — > 00 at the 
end of the calculation. 
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Let's discuss the properties of the operator H . Since \<f>\ — > 4>o as r ~ > °°j the continuum spectrum of H 1 - starts 
at eigenvalues |A| = m = \h\<pQ. H may also have bound states. We let mb be the smallest positive eigenvalue of 
H . By dimensional reasons, m?, ~ m. Now, observe, 

{a 3 ,H ± } = (16) 

Thus, a 3 maps a properly normalized eigenstate |A) of H with eigenvalue A into a properly normalized eigenstate 
of H ± with eigenvalue —A. Moreover, since a 3 maps zero modes of H 1 - into zero modes of H , all the zero-modes of 
H 1 - can be classified by their eigenvalue under a 3 . Writing, X(x) = (u(x), v(x)), we note that the zero modes of H 1 - 
with a 3 = 1 satisfy v = 0, T>'u = 0, while the zero modes of H with a 3 = — 1 satisfy u = 0, T>v — 0. So letting N+ 
be the number of a 3 = 1 zero modes, and iV_ the number of a 3 = —1 zero modes, we have, 

N = N+ - N_ = dim(ker(V i <)) - dim{ker(V)) = Index(H ± ) (17) 

Hence, N is the index of an elliptic operator, which is usually a strongly topological quantity. N has been first 
computed explicitly for a particular background string configuration in [13j to be: 

N = n (18) 

This result was later generalized 0i@ to arbitrary background string fields. 

We now return to the operator Hj~. Observe, [HkjH 1 - ) — 0. So, we can obtain the spectrum of Hk from the 
spectrum of H in the following way. Let, 

H x X(x) =AA(x) (19) 

First, suppose, A > 0. Then, the state, 

ip(x) = ciA(x) + c 2 a 3 \{x) (20) 
is going to be an eigenstate of Hj. with eigenvalue E, provided that, 

The eigenvalues of the above equation are, 



E = i^A 2 + fc 2 (22) 



and the eigenvectors, 



ci \ 1 / ±sgn(k)((X 2 + k 2 Y ± A) * 

C2 ^± (2(A 2 + fc 2 )^)" V ((A 2 + fc 2 )^ T A)" 



(23) 



Thus, each eigenstate of H 1 - with positive eigenvalue, generates one positive energy and one negative energy eigenstate 
of H^. However, this correspondence has to be taken with a grain of salt, since most eigenstates of H 1 - are continuum 
states, and the "1 to 2" map discussed above between eigenstates of H and eigenstates of need not preserve the 
density of states. 

The zero modes of H 1 - are also simultaneously eigenstates of Hk- These have the dispersion, 

E = ka 3 (24) 

So the zero modes of H become chiral fermions moving up or down the string depending on the sign of their 
eigenvalue under a 3 . 

C. Current - Naive Approach 

We now proceed to the computation of electric current at finite /i, T. This is given by: 

J 3 = e [ d 2 x($ 1 3 ^) = e I d 2 xtr{x\a 3 n(H)\x) (25) 
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where, 

n{E) = S 9< E ) (26) 
is the usual Fermi-Dirac distribution. Summing over each momentum sector fc, we obtain 3 , 

j3 = e lJ2 I d 2 xtr(x\a 3 n(H k )\x) = e~ Tr{c?n{H k )) (27) 

k ~* k 

Using the correspondence between spectra of H k and H^, we may schematically write the operator trace l|27|) as: 
j3 = e lY, E (iJE\a^ E )n{E) = e^Y.\ E (^,k, s \a 3 \^ s )n(E s (X,k))+ £ (A|a 3 | A)n(£(A, fc)) 

fe -E(-f/ fc ) fc \A(ff-L)>0,s=± A(//-L)=0 y 

(28) 

Here E(Hk) denote eigenstates of Hk, X(H ) denote eigenstates of -ff -1 , and ^A,fc,± denote eigenstates of generated 
by an eigenstate |A) of H- 1 , with energies E±(X,k) — ±\/ A 2 + A: 2 . Again, we stress that the above representation 
would have been absolutely correct if all the states contributing to the operator trace were discreet, and normalizable 
(for instance if T = and (i < m). In our case, this is not generally so, but we choose for now to ignore this problem, 
in order to illustrate the general idea behind the computation. We will later return to take the continuum states into 
consideration more carefully. 

For the moment suppose, T = 0, < /i < m&. Then n(E) — 8(E)6(fi — E). Hence, only states generated by zero 
modes contribute to the sum in J5HJ), as all the other states have energies \E\ — \/A 2 + fc 2 > |A| > m& > /i. The zero 
modes are eigenstates of a 3 , and thus, satisfy, E — a 3 k and (A|a 3 |A) = a 3 . Thus, 

j3 = 4 E^+^ " fc W fc ) - N - 6 ^ + k W~ k » = eN I f- 9 ( k )°^ -k) = ^ = ^p (29) 

k J IT IT IT 

where we've used the fact that index TV is equal to the winding number of the vortex n. 

Now, let's relax our assumption and work at arbitrary T, /i. We first need to evaluate the matrix element 
(V'A,fc,s|a 3 |V'A,fc,s) (in what follows we suppress the indices A, k,s). Using ea. i(2U|) and (a 3 ) 2 = 1, we see (%p\a 3 \ip) = 
(|ci| 2 + |c 2 | 2 )(A|a 3 |A) + (c*c 2 + c* 2 c x ). Recalling a 3 |A) = | — A), we obtain, {ip\a 3 \ip) = c\c 2 + c* 2 ci = -|, where we've 
used eq. Hence, 

j3 = e |E( E Y^yME^k^+N+nW-N-ni-k)) (30) 

k \X(H ± )>0,s=± ' ' J 

Now, observe that E S (X, k) = E S (X, —k) for A > 0. Hence, the first sum in the brackets in eq. (|30|l is odd in fc, and, 
thus, the contribution to J 3 from non-zero modes of H 1 - cancels out exactly, leading to: 



J 3 = eN— ^ n ( k ) = en [ —n(k) = 

7„ J 



enn (fj,,T) (31) 



Here no(/i,T) is the number density of a free massless chiral fermion in 1 dimension, at finite chemical potential /i 
and temperature T. It is a peculiar fact that no(n,T) is temperature independent and equals ■Jjr-, so that, 

= = (32, 



D. Current - Corrections from Polarized Continuum 



Although, the result i|32|) is very attractive, it is actually generally incorrect. We know that this result is exact for 
T = 0, jj, < to, when J receives contributions only from normalizable eigenstates of Hk- We will now show, that 



Here tr denotes matrix trace and Tr denotes a general operator trace. 
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the presence of long range vortex fields polarizes the continuum eigenstates of H k in a way, which might significantly 
modify the result l|32[l for /i > m. 

Let's return to the trace (|27|l . We can rewrite this expression in terms of spectral current density as: 

J 3 = J dEn(E)f(E) (33) 

f(E) = ei-^TTrV^-S)) (34) 

fc 

We use the following representation of the delta function, 8{x) — ^ lim e ^ + (:jqbe — J^u) > ^° rewr ite, 

f{E) = A lim e\YTr ( a 3 ( TT 1 - - — )) (35) 

J v 7 2tt £ ^o+ \ y H k + z+ H k + z-'J y ' 

where z + = —E + ie, z~ = —E — ie. From here on, the limit e — > + is implied. Simplifying Ij35(l . 

■3f F \ _ I'^if.T f 3/ #fc^ + H k -z~ \ _ lr^ie / 3 ka 3 + H ± -z+ kg 3 + - z~ ' 

1 j L^ 27r V l ff 2 -(z+) 2 Hl-{z-f ] ) L^2n r { a V 2 + fc 2 - (z+) 2 H^+k 2 -(z-) 



2 ' 



(36) 

where we've used = if + fc 2 . The terms in <|36(1 . which are odd in fc cancel out, and we obtain, 



' 3 (.El ^ * 6 T \ 3 ( ~ Z+ ^ ~ ' Z m 

)_ L^27 r ^ ( i/^ 2 + fc2 _ ( z + )2 ~ #± 2 - fc2 - ( z -)2^ 

Now, {a 3 , ff- 1 } = 0. Hence, for any function /, Tr{a 3 H jL f{H 1 - 2 )) = -Tr{H^ f(H J - 2 )a 3 ) = -Tr(a 3 H x f(H^ 2 )) = 0, 
and, 

f {E) = IV ^Tr f a 3 ( ^— — ^— ^ ) ] (38) 

L fc 2n V ^ +fc 2 -(z+) 2 H ±2 + k 2 -(z-) 2 I 

We now introduce the function g, 



^ = Tr (a 3 f ) = Tr ( J] - Tr ( *** \ (39) 

This function is very well known 0j@ as it satisfies, 

N = Index{H^ ) = Jim ^g(M 2 ) (40) 

More generally, g(M 2 ) is related to the spectral asymmetry a k (E), of the Hamiltonian H k 7], and hence to its 
^-invariant as, 

a k (E) = l-k{G{k 2 -{z+) 2 )-G{k 2 -{z-) 2 )) (41) 

r, k = 2 / dEa k {E) (42) 
Jo 

G(z) = ^ (43) 

Here, g(z) is understood as the analytic continuation of g from 5R+ to C. From eq. H38JI . we can express j 3 (E) in 
terms of G as, 

?w = z E ( Z+G(fc2 ~ (z+)2) - z ~ G{k2 ~ ( ° 2) ) = I £ e f afe(i?) (44) 
& fc 
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Following the technique of trace identities described in detail by one can explicitly calculate g(M 2 ) to be: 

M 2 

g{M 2 ) = n~{n-<S>) (45) 
777 + M z 

Hence, the index N — lim M 2^ g(M 2 ) — n, in agreement with previous calculations Continuing g to the 

complex plane, we obtain, 

(*) = 2_ (n _*)_L_ (46) 

Hence, generically, G has a pole at z — 0, and a pole at z — — m 2 , i.e. at the continuum threshold. Notice, however, 
that the pole at z = m 2 disappears when n = 

We can now substitute the result (|46|l into l|41|l and take the limit e — > + to calculate the spectral asymmetry, 

a k (E) = k sgn(E)(n 6(E 2 - k 2 ) - (n - <P)5(E 2 - k 2 - m 2 )) (47) 

which yields the 77- invariant, 

k 

r) k =n sgn(k) - (n - $) (48) 

We note that eq. I|48|l is in agreement with previous calculation of the 77-invariantQ. 

Returning to the evaluation of current, we substitute the result 147|) into eq. 144fl to obtain, 



f 



{E) = ~ I (n(6(E -k) + S{E + ft)) - (77 - - V^ 2 + ™ 2 ) + S(E + yjk 2 + m 2 ))) (49) 



and the total current in the string direction becomes, 



2vr 

This can be conveniently rewritten as, 



J :! < I * ^77 77(fc) - "^—^(ni^k 2 + m 2 ) + n(-^k 2 + m 2 ))^ (50) 



J 3 = e(nn (n,T) - - - ' n m (n,T)) (51) 



where no(fi,T) = £- is the familiar number density of one-dimensional chiral massless fermions, and, 

— (n(Vk 2 + m 2 ) + n{-^k 2 + m 2 )) (52) 

27T 

is the number density of one-dimensional 2-component (Dirac) fermions of mass 777. 

Several comments are in order here. First of all, we see from eq. H51(l that the naive result 1)31(1 is generally modified 
by a contribution from modes located at the continuum threshold. Observe, that for ji = 0, the current J 3 vanishes 
for all temperatures. At non-zero chemical potential, two cases are of particular interest. The first case is that of a 
local string, satisfying the finite energy condition, Dtp — * faster than 1/r, which implies <& = n. In this case, the 
contribution from continuum modes vanishes, and we recover our initial result Il.'j2[| . which is due solely to the zero 
modes, 

This "coincidence" can be explained as follows. If Dip — > fast enough, the fields in the problem are, essentially, 
short range, and hence we can easily put the system in a box, making the spectrum discrete, so that the argument in 
section II C is correct. 

Let's briefly discuss what happens when we add the second fermion ip to the problem. Recall, we used this fermion 
to cancel gauge anomalies of our model. As noted in section II A, the contribution of to J 3 can be obtained by 
taking e — > e, /1 — > /t, 77 — > —77, $ — ► —<fr. In particular, the continuum modes at threshold again cancel out, and 
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J 3 = — In particular, if the chemical potentials of ip and ip fermions are the same, we can obtain a non- vanishing 
total electromagnetic current along the string, by letting 4 e = — e , so that, 

4m = ^ (54) 

The second practically interesting case is that of a global string. This case can be recovered by taking $ — * 0. 
Then, 

J 3 = en(n (M, T) - in ro (/x, T)) (55) 

In this case, the field <fi is long range, and there is a significant modification of the result J3U- Note that n m (/j,,T) is 
no longer temperature independent, so for simplicity, we choose to work at T — 0, \i > 0. Then, n(E) = 9(E)d(fi — E) 
and, 



J 3 = — (u - 0{u - m) Vm 2 - ™ 2 ) (56) 

Z7T 

Thus, for /i < m, the current is governed by our original result (|32|l . while for /i > m, we also get a counterflow current 
from the states at continuum threshold. Thus, J 3 (^) has a cusp at fi = m, and for fi ^> m falls off to as ~. 



III. CONCLUSION 



In this paper, we have found an exact expression for the electric current on superconducting strings as a function 
of fermion chemical potential and temperature. We've analyzed the case of both local and global strings, and our 
analysis has not been limited to a low energy theory of zero modes in the string core. Our ability to obtain such an 
exact result has been due to a cancellation (or partial cancellation) of contributions of all, but the zero fermion modes 
to the current. For local strings, we've seen that for all values of T and /j, the current is due to zero modes in the 
string core. On the other hand, for global strings, the current receives contributions both from the zero modes and 
from certain states at continuum threshold. The latter contribution tends to cancels out the contribution from the 
zero modes as the fermion chemical potential becomes much larger than the fermion mass m. The results for /i»ra, 
might be particularly interesting in application to currents on axial strings in dense quark matter [T^. where the gap 
A < fi. 

We would like to note that the study of persistent topological currents and spin currents in conjunction with 
problems, such as, for example, Quantum Hall Effect [T^| and Spin-Hall Effect [Tlj. has over the past years become 
an active subject of research in condensed matter physics. It would be interesting to investigate the relation of the 
phenomenon discussed in this paper to problems in condensed matter systems. For instance, persistent supercurrents, 
are known to appear on vortices in superfluid 3 He — A and, somewhat similarly to currents considered in this paper, 
are due to chiral anomalies 16J. 
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